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1 Introduction 



In this paper we find some self-similar exact solutions of the following version of the 
Laplacian growth problem in the half plane H introduced in pQ. Let 7 be a smooth 
non-self-intersecting curve in HI from a point x_ G R to a point x + G R (we assume that 
x_ < x + ) and B be the domain bounded by this curve and the segment of the 

real axis. (In [1] such domains are called fat slits.) Suppose the curve 7 = 7(T) moves 
with time T according to the Darcy law: 

V n (z) = ~d n <f>{z), 267. (1) 

Here d n is the normal derivative at the boundary, with the outward looking normal vector, 
V n (z) is the normal velocity at the point z G 7 and is a unique harmonic function in 
H \ B such that 

(i) = on 7 and on the rays of the real axis [—00, [x + , +00]; 

(ii) <p(z) = Im z + o(l) as Im z — > +00. 

As is argued in [T], the growth process is well defined if both angles a± between 7 and 
the real axis at the points x± are acute. Then these angles as well as the points x_, x + 
stay fixed all the time. 

Comparing this setting with the standard Laplacian growth in the upper half plane 
(see, e.g., [2]), we see that the conditions on the harmonic function are very similar: 
= on an infinite contour from left to right infinity and becomes Im z as Im z — > +00. 
An important difference is that in our case, unlike in the standard one, only a finite part 
of the level line = (namely, the part which lies above the real axis) moves according 
to the Darcy law while the remaining part (the rays of the real axis) is kept fixed despite 
the fact that the gradient of is nonzero there. 

It is often convenient to treat the growing domain B as an upper half of the domain 
D = B U B symmetric with respect to the real axis, with the boundary T = 7 U7. Here B 
is the domain in the lower half plane obtained from B by complex conjugation z —>■ z. In 
what follows we call such domains simply symmetric. Then one can extend the problem 
(CE]) to the whole plain "by reflection", i.e., by saying that V n (z) = \ d n 4>(z) for z G 7 
and V n (z) = — | d n (p(z) for 267, where is a harmonic function in C \ D such that (i) 
(j)(z) = — 4>{z); (ii) = Oonr = 7U7; (iii) <j>{z) = Xmz + o(l) as Xmz — > ±00. 

2 Formulation in terms of conformal maps 

As is customary in moving boundary problems, we reformulate the problem in terms of 
time-dependent conformal maps to (or from) some fixed reference domain from (or to) 
the domain C \ D(T), where the Laplace equation is to be solved. In our case there are 
two distinguished choices of the reference domain: the upper half plane and the exterior 
of the unit circle. 
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Figure 1: The conformal map p(z). The mirror image of the domain B in the lower 
half-plane is also shown. 

The upper half plane. This choice is natural when one deals with the original for- 
mulation of the problem in the upper half plane not extending it to the lower one by 
reflection. Let p(z) be a conformal map from H \ B (in the "physical" z-plane) onto H 
(in the "mathematical" p-plane) shown schematically in Fig. [TJ We normalize it by the 
condition that the expansion of p(z) in a Laurent series at infinity is of the form 

p(z) = z + — + u k z~ k , \z\-kx (2) 

Z k>2 

(a "hydrodynamic" normalization). Assuming this normalization, the map is unique. 
The upper part of the boundary, 7, is mapped to a segment of the real axis [p_, p + ], 
while the rays of the real axis outside B are mapped to the real rays [—00, p_] and 
[p+, 00]. From this it follows that the coefficients u k are all real numbers. The first 
coefficient, u%, is called a capacity of B. It is known to be positive. We also need the 
inverse map, z(p), which can be expanded into the inverse Laurent series 

z(p) =p + akP~ k , \p\ -> 00 (3) 

with real coefficients a k connected with u k by polynomial relations. The series converges 
for large enough \p\. 

Clearly, the harmonic function Imp(z) obeys all the conditions required from the 
<p(z), so (j)(z) = Imp(z). It is easy to see that d n Xmp(z) = \p'(z) \ on 7, so one can write 
the Darcy law as V n (z) = hp'(z)\. 

The Schwarz symmetry principle allows one to extend this reformulation to the whole 
plane as follows. The function z(p) admits an analytic continuation to the lower half 
plane as zip). The analytically continued function performs a conformal map from the 
whole complex plane with a cut on the real axis between p_ and p + onto the exterior 
of the symmetric domain D = B U B. Correspondingly, the inverse function, p(z), obeys 
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p(z) = p(z) and d n Tmp{z) = — \p'(z)\ on 7. The evolution of the whole closed contour 
T = 7 U 7 can then be written in a unified way as V n (z) = ^\p'(z)\ for any z G T. 

The map z(p) plays a crucial role in the embedding of the problem into the disper- 
sionless KP hierarchy found in [1] but appears to be rather inconvenient for constructing 
explicit solutions. This is certainly related to the fact that the symmetrically extended 
reference domain in the mathematical plane is singular (plane with a cut) and, moreover, 
depends on time. There is another choice of reference domain which seems to be less 
natural from the point of view of integrable hierarchies but is more suitable for finding 
explicit solutions. 



The exterior of the unit circle. In this case it is convenient to work with sym- 
metrically extended domains from the very beginning. Let f(w) be the conformal map 
from the exterior of the unit circle onto the exterior of the symmetric domain D such 
that /(oo) = 00 and f'(oo) = r > 0, so that the Laurent expansion at 00 has the form 
f(w) = rw + uq + 0{w~ 1 ) with real coefficients. The coefficient r is called the (external) 
conformal radius of the domain D. The connection with the previously considered map 
is as follows: let w(z) be the inverse function, then p(z) = r(w(z) + w~ 1 (z)) + u with 
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\(P+ ~ P-) and «o = \(j>+ + P- 



Let us rewrite the Darcy law as a dynamical equation for the function f(w) = f(w, T). 
To do that, we need the following simple kinematical relation which can be derived in a 
direct way. Let z(a, T) = x(a, T) + iy(a, T) be any parametrization of the contour such 
that a is a steadily increasing function of the arc length, then 

V n = ^{d T xd a y - d T yd a x) = ^-(d T zd a z - d T zd a z) , (4) 
al Aim 



where dl = \dz\ = y (dx) 2 + (dy) 2 is the line element along the contour. According to our 
convention, V n (z) is positive when the contour, in a neighborhood of the point z, moves 
to the right of the increasing a direction. Set o = —ilogw, z = f(w,T), z = f(l/w,T), 
then the kinematical formula gives 

Vn(f(w)) = — ^ (d T f(l/w)d w f(w)-d T f(w)d w f(l/w)) , \W\ = 1 

On the other hand, the Darcy law reads 

, r , ,/ ^ 1 \dp/dw\ 

K(/W) = 2WAM' H 

Equating the right hand sides and using the relation dp/dw = r(l — w~ 2 ), we get the 
equation 

wd w f(w)d T f{l/w) - wd w f(l/w)d T f{w) = r w-w' 1 , \w\ = 1 (5) 

which is the main dynamical equation of the problem in terms of the conformal map. 

It is instructive to compare it with the similar equation for the Laplacian growth 
process in the whole plane (i.e., without the condition that <fi vanishes on the the real 
axis) and with the same type of source at infinity. In the latter case the evolution 
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necessarily destroys the reflection symmetry, so the function f{l/w) should be replaced 
by f(l/w) = f(l/w) while the form of the r.h.s. should be also changed to — ir(w — w^ 1 ). 

For finding explicit solutions of prime importance is the case when the function 
f(w,T) admits an analytic continuation across the unit circle for all T in some time 
interval, so that f(w,T) is actually analytic not only in its exterior but in some larger 
domain containing it. However, one should take into account that such a continuation 
is impossible through the points w = ±1 which are pre-images of the two corner points 
and so f{w) is not analytic there. Assuming that it is analytic everywhere else on the 
unit circle, one can analytically continue equation ([5]) as follows: 

!— ir{w — w^ 1 ) , Imw > 
(6) 
ir{w — w^ 1 ) , Imw < 

Note that the analytic continuations to the upper and lower half planes are different. 



3 Formulation in terms of Schwarz function 

It is known [2J that the standard Laplacian growth problem in the whole plane can be 
integrated in terms of the Schwarz function. For analytic contours, the Schwarz function 
S(z) is defined as the analytic continuation of the function z away from the contour. In 
other words, the Schwarz function for a curve F is an analytic function S(z) such that 
z = S(z) for z G T (see [3] for details). If the curve depends on time, so does its Schwarz 
function, S = S(z,T). 

In our problem the contour V is not analytic because it always contains the two corner 
points. However, if the curve 7 is analytic, its Schwarz function, which we still denote 
by S(z), is well defined in a strip-like neighborhood of the curve, with the width of the 
strip tending to zero around the endpoints x±. We thus have z = S(z) for z G 7 or, 
equivalently, z = S(z) for z G 7- The latter formula just means that S(z) = S(z) is the 
Schwarz function for the complex conjugate curve 7- Therefore, we can expect that our 
problem can be integrated in terms of two Schwarz functions, S(z) and S(z), enjoying 
equal rights. Below, just for brevity, the pair (S(z), S(z)) is referred to as the Schwarz 
function of the piecewise analytic contour T = 7 U 7- In the mathematical w-plane, the 
role of the Schwarz function is to connect f{w) and f(l/w) in their common domains of 
analyticity: 

f S(f(l/w)), lmw>0 

fW = (7) 
[ S(f{l/w)), lmw<0 

Note that the second equality (at Imw < 0) is obtained from the first one by complex 
conjugation. 

In complete analogy with the Laplacian growth problem in the whole plane, the 
growth process flTJ is encoded in the equation id^S = d z p in the upper half plane. 
Extending it to the lower half plane by reflection (i.e., complex conjugation), we have, 
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for all^eC\D(T): 

' id T S(z,T) = d z p(z,T) , lmz>0 

(8) 

t id T S(z, T) = -d z p(z, T), lmz<0 

It is convenient to rewrite the above relations in terms of the piecewise analytic 
function 

f i(S(z) - z), lmz>0 

M(z) = (9) 

[ -i(S(z) -z), lmz<0 

introduced in pp. By construction, it provides the analytic continuation of the function 
2\1mz\ away from the contour V. Equivalently, we write S(z) = z — iM(z) for Imz > 
and S(z) = z + iM(z) for Im z < or, passing to the mathematical w-plane, 

( f(l/w) + iM(f(l/w)), lmw>0 
f(w) = (10) 
[ f(l/w)-iM(f(l/w)), lmw<0 

In terms of the function M = M(z, T) the dynamical equation ([8]) acquires the form 

d T M(z,T) = d zP (z,T) (11) 

Note that the function d z p(z, T) is holomorphic in C\D(T) with the expansion d z p(z, T) = 
1 + 0(1/ z 2 ) at infinity 

The function M(z) can be uniquely decomposed as M(z) = M + (z) — M_(z), where 
M + is analytic in D and M_ is analytic in C \ D with zero at infinity. This decomposition 
is given by the integral of Cauchy type 

.,- / n 1 / M(z')dz' 1 / \lmz'\dz' ^ , 

MAz) = & —-^ = — ? 1 , zeD 12 

2ni JdD Z' — Z Til JdD z' — z 

(and the same integral for M_(z) with z outside D). The function M + is analytic every- 
where in D and, by our assumption, can be analytically continued across the arcs 7 and 
7 everywhere except for their endpoints on the real axis, where it has a singularity. One 
can see that the function M + (z) is the generating function of integrals of motion for our 
problem. Indeed, it is straightforward to calculate its time derivative: 

dT M + (z) = ^<f ggg^fWgO Wl (13) 

7TZ JdD Z — Z' 

(see [1] for details). Plugging here V n (z) = i|p'(2;)| and recalling that dp(z) is purely real 
on T, we obtain 

where only the residue at infinity contributes because p\z) is holomorphic in the exterior 
of D and z is inside. Therefore, we have obtained the important equation 

d T M + (z } T) = l (14) 
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which allows one to construct an infinite series of integrals of motion by expanding M + (z) 
into a series in z. 

In the case of general position, when X- ^ x + , one may expand M + (z) in a Taylor 
series around a point on the real axis lying inside D on the segment [x_,a; + ]. Choosing 
the coordinate in such a way that x_ < < x + , one can expand around the origin: 

oo 

M+(z) =T + Y,kT k z k - 1 . (15) 

k=2 

Equation (jXlj) shows that the coefficients T k , k > 2, are conserved: drTf. = 0. They are 
the following harmonic moments of the domain H \ B: 

T k = 4r i \Xmz\z~ k dz = — [ Xm (z~ k ) d 2 z, k>2 (16) 

Tlik JD 71K JH\B V ' 

One may also consider degenerate configurations with x- = x + = 0. In this case the 
function M + (z) — T is still conserved but can not be represented by a Taylor series around 
the origin, and so the moments ffTBT) are ill-defined. 



4 Self-similar solutions 

Self-similar solutions are characterized by the property 

f(w,T)=r(T)f(w) (17) 

which means that the time evolution is equivalent to a dilatation in the z-plane. Because 
the time evolution is always such that the points x± stay fixed, self-similar solutions have 
chance to exist only when x_ = x + = which is just the degenerate case mentioned at 
the end of the previous section. Exact self-similar solutions of Laplacian growth in radial 
and wedge geometries were studied in 



4.1 General relations 

Differential equation for f(w). Substituting the self-similar ansatz (1X7)1 into (jUJ), we 
see that the variables separate: r(T) —T/A and 



—i[w — w 



_1 ) , Xm w > 



wd w f(w)f(l/w) -wd w f(l/w)f(w) = A 



i(w — w 



_1 ) , Xm w < 



where A is a constant. This latter equation can be treated as the Wronskian relation for 
a second order differential equation of the form 

A 2 (w)f + A^wW + A (w)ip = 

Namely, let ip2 be two linearly independent solutions, then their Wronsky determi- 
nant, W = ipi4>2 ~ ^2^1 obeys the relation 

d w \ogW{w) = - 



A 2 (w) 



7 



© 








Figure 2: A corner point at the origin. 



In our case W = ±iA(l 



w 



), so d w \ogW 



1 1 2 . , 

H (note that this ex- 

1 w 



w + 1 w 

pression is the same in the both half planes) and the functions f(w) and f(l/w) are two 
solutions of the ordinary differential equation 



d 2 f- 



1 



1 



W 



w 



with yet unknown function V(w), or 



w 2 dlf 



2w 



XV, 



-)d w f + V(w)f 
w J 



d w f + V(w)f = 







(19) 



(here V(w) = w 2 V(w)). 

From the requirement that f(l/w) is the second solution of equation ([191) it follows 
that V{w^ 1 ) = V(w) and from the behavior at infinity f(w) = w + 0(1) that V(w) = 
0(w~ 2 ). More precisely, if f(w) = w + ao + aiw' 1 + 0(w~ 2 ), then a simple calculation 
shows that V(w) = 2(1 — ai)/w 2 + 0(w~ 3 ). Further, from the fact that f(w) is regular 
outside the unit circle (except for the simple pole at infinity) and V^w^ 1 ) = V{w) it also 
follows that V(w) may have singularities only on the unit circle. The argument is the 
same as for self-similar Laplacian growth in the wedge [6J. Suppose, for example, that V 
has a pole at some point £ outside the unit circle. Then the only way to compensate it is 
to impose the condition /(£) = which means that the point £ actually lies on the unit 
circle because the only possibility for the point to be the image of any point w such 
that | w\ > 1 is to be a singular boundary point in which case its pre- image must belong 
to the unit circle. 

Here we do not address the question what is the most general type of singularities 
allowed for the function V(w) on the unit circle and will restrict ourselves by the case 
when V(w) is a rational function. Suppose the boundary has a corner with an angle 
9 when w moving along the unit circle passes through a point £ ^ ±1 (Fig. [2]), then 
in a vicinity of £ the conformal map has the form f(w) oc (w — £) 6,/ ' 7r . (We assume 
that the boundary can be non-smooth only when it crosses the origin, then /(£) = 0.) 
Substituting it to equation (fl9|) . we get 



- i--k 2 , 

(w - 4) v 



w 



(20) 



s 



Note that the main part is the same for the angles 6 and it — 6. Clearly, a similar pole with 
the same 6 is at the complex conjugate point £. Poles at £ = ±1 should be considered 
separately. In this case we get 

V(w)= y Tl) V (l + OfoTl)), ™-±l (21) 

where # is now the angle between the real axis and the nearest piece of the boundary 
curve emanating from the origin. 

Let us summarize the properties of the function V(w): 

a) V(w) a rational function of w such that V(w) = V(w~ 1 )] 

b) The only singularities of V(w) are second order poles on the unit circle with the 
principal part given by (120]) . ( l21i) : 

c) V(w) = 0(w~ 2 ) as w — > oo. 

In some cases these conditions allow one to fix V(w) uniquely (see the examples below). 
In order to find the conformal map f(w) one should solve equation ffl9|) . choose a solution 
such that f(w) = w + 0(1) as w — > oo and check that the map is indeed conformal, i.e., 
that /'(n;) 7^ for all w outside the unit circle. 

In fact, it is enough to check that f(e tip ) makes just one complete turn along the 
boundary of D in the counterclockwise direction as (p runs from to 2ir (the argument 
principle, see, e.g., [H)])- Indeed, the number of zeros of the analytic function f'(w) 
outside the unit circle is given by the integral — ^ §i w i =1 daxgf'(w). (More precisely, to 
regularize possible singularities in pre-images of the corner points one should consider 
the integral — ^ §\ w \= e * dargf'(w) with a positive e — > 0.) On the other hand, since / 
is conformal in a vicinity of the unit circle, the normal vector to the boundary swings 
through the angle arg /'(«;) under the map, i.e., arg/'(w) = 4>(f(w)) — axgw, where 
is the angle between the outward pointing normal vector to the boundary of D and the 
real axis. From this it is clear that arg/'(w) does not change when w makes a round trip 
along the unit circle and thus f'(w) can not vanish outside it. 

The form of M + (z). For the conformal map p(z,T) self-similarity means p(z,T) = 
r(T)p(z/r(T)) and analogously for the M-function: M(z,T) = r(T)M(z/r(T)). Plug- 
ging this into ffTTj) . we get r = 1/A, i.e., r = T/A and 

M(z) - zM'(z) = Ap'(z) 

The function M(z) here has the meaning of the M-function for the domain with conformal 
radius r — 1, i.e., at T = A. Integrating both sides along the boundary with the Cauchy 
kernel, we get the following condition for M + : 

M+(z) - zM' + (z) = A 

with general solution M + (z) = cz + A, or 

M+(z,T) = cz + T (22) 

where c is an arbitrary real number. It is important to note that in the case under 
discussion the domain D consists of at least two disconnected pieces whose boundaries 
intersect just at the origin, and the constant c can be different in the different pieces. 
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4.2 Examples 

One-petal solution. The simplest example is the one-petal pattern growing in the 
upper half-plane and symmetric w.r.t. the imaginary axis (Fig. [3]). Let a be the angle 
between the tangential lines at the origin and the real axis, then the conditions a), b), c) 
above fix the function V(w) uniquely: 

Tr , , 16a / a\ w 2 . „, 

V(w) = 23 

71 V 7T J [W 2 — 1) Z 

After the substitutions w = e x , f(w) = (smhx) 1 ^ 2 ip(x) equation ffl9|) takes the form 



r+^rj^-^-U (24) 

sinn x 4 



where 

2a 1 1 1 

7 = V"2' -2 <7< 2' (25) 
Its solution having the required properties at infinity is 



^(x) = e x / 2 (l-e- 2x y 2 F l (^ in i 



7,7-1/2 
/2 



e~ 2x 



The hypergeometric function with these parameters can be expressed through elementary 
functions: 

2F 1 ( 7 ' 7 /2 1/2 ; z 2 ^ = \ [(1 + z)^ + (1 - (26) 
For /(w) we thus obtain: 



/H = -Vw 2 - 1 [(1 - w- 1 ) 7 (l + w- 1 ) 1 " 7 + (1 + ^-^^(l - w- 1 ) 1 -^ (27) 



At 7 = (a = 7r/4) this gives the solution /(w) = \Jw 2 — 1 mentioned in [T]. It is the 
Bernoulli lemniscate given by 

(x 2 + 2 / 2 ) 2 = 2( 2/ 2 -a; 2 ) (28) 

in Cartesian coordinates or r 2 = —2 cos(26 l ) in polar ones (Fig. HI dotted line). Using the 
argument principle, one can show that f'(w) ^ at \w\ > 1 for all 7 between —1/2 and 
1/2, i.e., we obtain a continuous family of one-petal self-similar solutions for all angles 
< a < 7r/2. The examples for a = tt/8 and a = 3tt/8 are shown in Fig. |H 



One-petal solution via integral equation. The same results can be obtained in a 
different way, using an integral equation for the conformal map f{w) which follows from 
the properties of the M-function. The function M + for the one-petal solution has only 
one cut which is the whole real axis and by virtue of fT22|) the jump across the cut must be 
a linear function. The coefficient can be fixed by analyzing the local behavior around the 
origin. Indeed, let e ie ^ = dz/\dz\ be the unit tangent vector to the curve (represented as 
complex number) at a point z, then, by definition of the Schwarz function, S'(z) = e.- 2lQ i z ) 
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(a) 



(b) 



Figure 3: (a) One-petal solution, (b) two-petal solution. The mirror images in the lower 
half-plane is also shown. 



for z on the curve in the upper half plane and S'(z) = e~ 2%e ^ for z on the curve in the 
lower half plane. These formulas allow us to calculate discontinuity of the function M' 
(which obviously equals the discontinuity of M' + ) across the real axis. Using (jUJ) we find: 

[M'\ := M' up - M' down = h2-S' up - S' down ) = ~(2 - e~ 2ia - e 2ia ) = -4z sin 2 a (29) 

i V 

Thus in the upper half-plane 

M+ = -2i sin 2 az + A, (30) 

and in the lower one 

M+ = 2i sin 2 az + A, (31) 

Using ([7]) we can express the sum of the values of f{w) on the opposite sides of the 
cut from —1 to 1 in the w-plane. Set {f{w)} = f(w + iO) + f{w — iO), then 

{f(w)} = 2f(l/w) + i[M + {f(l/w))} = (2 - 4 sin 2 a)f(l/w) = 2 cos 2a f(l/w) 

where it is taken into account that 

[M+(f(w))} = M + {f(w + iO)) - M+{f(w - z0)) = -4i sin 2 af(w). 

We have obtained the equation 

{f(w)} = 2cos2af(l/w), w G [-1, 1] (32) 

which determines the conformal map. 

It can be transformed to an integral equation by means of the substitution 



#0) 



' Vw 2 - 1 
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-1.0 




1.0 



Figure 4: One-petal solutions: a = | (dashed line), 
a = | (dotted line) and a = ^ (solid line). 



which allows one to transform the mean value of the function / on the cut to the jump 
of the function g. Note that the branch of the square root should be chosen such that 
= wy/1 — w~ 2 be an odd function outside the unit circle. The symmetry of the 



w- 



petal implies that f(—w) = —f(w), thus g( 
equation (1521 acquires the form 



-w) 



g(w). In terms of the function g 



[<?H] 



-2i cos 2a 



g(l/w) 



w 



W E [-1, 1] 



(33) 



The function g is analytic everywhere in the w-plane except for the two branch points at 
w = ±1, with the jump across the cut [—1, 1] being given by (1331 . Taking into account 
that g(oo) = 1, we can represent it by an integral of the Cauchy type along the cut: 



9{w) 



1 



2i cos 2a 



g(l/x)dx 



2m J -i x(x — w) 

Finally, using the fact that g is an even function, we arrive at the integral equation 

2cos2a f 1 g(l/x)dx 



g{w) 



TV 



x 2 — w 2 



(34) 



We know that f(w) oc (w — 1) " in a vicinity of w — 1, so we need a solution such that 

2a. 1 

g[w) oc (w — 1)~~2 near w — 1, and similarly for a vicinity of w — — 1. Note that the 
exponent is the same 7 as in (|25j) . Consider first the case of positive 7. After a further 



substitution g(w) 



w 



-2\ 7 



G(w ), using the fact that #(±1) = (valid for positive 



7), we arrive at the equation 

(1 -xy- l G(x) 



sin7T7 r 1 r ! / 2 (l - ty~ x G{t)db 

TT Jo 1 — Xt 



(35) 
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Comparing it with equation (42) from |9j, we can immediately write down the solution: 

G(i)=^( 7 'V /2 ;*), 

which is the same function as in ( |26l) . The map f(w) is then given by the explicit 
formula ( 1271) . A direct substitution shows that the function g(w) = f(w)/\ / w 2 — 1 obeys 
the integral equation (134]) for both positive and negative 7, I7I < 1/2. 



Two-petal solution. Consider now a two-petal pattern growing in the upper half- 
plane and symmetric w.r.t. the imaginary axis (Fig. [3]). Let the angle between the two 
petals be 2/3. Clearly, the angles satisfy the condition < (3 < n/2 — a. Now V{w) has 
4 double poles. Two of them are at w = ±1. The symmetry implies that the other two 
poles are at w = ±2. The function V(w) is again uniquely determined by the conditions 
a), b), c): 



V(w) 



16a 



7T 



1-2 



w 



2/T 



— 1- — 



w 



7T 



7T 



(w 2 + iy 



(36) 



nj (w 2 -l) 2 

Note that the same function corresponds to the angle ir/2 — (3. The same substitution 
as above, w = e x , f(w) = (sinh x) 1 ' ' 2 'ip(x), brings equation ( 1T91) to the form 



„ | 7(1 - 7) ^ gl ~ 5) , _ 1 , 
sinh 2 x cosh 2 a; 4 

where 7 is as in (125!) and 5 = 2(3 /it. A further change of variables, 



if) = 2 (sinh x) 7 (cosh x) 2 7 F(l/cosh 2 x) , t 



1 



cosh x 



puts the equation in the canonical hypergeometric form 

t(l - t)F tt + + j)t\ F t - l((2 T +25-l)(2 T -25+l)) J P = 0. 
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Two linear independent solutions are 

(2 7 +2<J-l)/4, (2 7 -25+l)/4 



F w (t) = 2 F X 



and 



F {2 \t)= t l '\F 1 



1/2 



/ (2 7 + 2<5+l)/4, (2 7 -25+3)/4 



V 



3/2 



For the function if> the first solution gives 

/ (2 7 +25-l)/4, (2 7 -25+l)/4 



ip(x) = (sinh x) 7 (cosh x) 



V 



1/2 



; cosh x 



(37) 



(3f 



13 



Passing to the function f(w) and to the original angles a, (3, we have: 



2\1- 



f(w) = w(l - w~ 2 )-(l + w~ 2 ) 



2 Fi 



71 



71 



V 



1/2 



(39) 



/ 



The symmetry implies that the conformal map must be an odd function of w. The 
function given by (1391) is indeed odd while the second solution, F^ 2 \ leads to an even 
function of w. Using the argument principle, one can see that for 0<(3<a<ir/2 the 
derivative f'(w) has no zeros in the exterior of the unit circle. Therefore, f(w) given by 
f l39|) is the conformal map for the two-petal solution. 

Note that this function looks somewhat simpler in the variable p = w + w" 1 living in 
the upper half plane: 



r(p) = p(l - 4p~ 2 ) a ^ 2 Fx 



( a+(3 _ a-P 



71 



71 



4p~ 



\ 



1/2 



/ 



(40) 



The analytic continuation of the function (HUT) to the region \p\ < 2 reads: 



z{p) = 2(l-p 2 /4) a/7T 



le 



-i/3 



r(^)r( 



' 7T 2' 7T 



7T " 1_ 2 



T 



+ e 



r(|)r. 



2/3 



P 



1-2/3/tt 



2^1 



' 7T ~ l ~ 2' 7T 



(41) 



P 



V 



3_ 2/3 

2 7T 



The boundary of the two-petal pattern in the upper half plane is obtained as z(p + iO) 
for -2 < p < 2. 

In the case (3 = a the hypergeometric function in fj4*0|) becomes trivial (equal to 1). 
The corresponding conformal map appears to be a map to a slit domain: the function 
z{p) = p(l — 4p~ 2 ) a / 7r takes the upper half plane to the upper half plane cut along two 
straight segments emanating from to the NE and NW quadrants at the angle a to 
the real axis. In fact this means that (3 = tt/2 — a rather than a but as it was already 
mentioned, the function fl36l) is the same in both cases. 



Some typical two-petal solutions are shown in Fig. [51 The following cases are to be 
distinguished: 



A) a e (0, 1) 

i) (3 G (0, a) - the petals are shown in Fig. [5] (a), they become very thin as (3 
approaches a from below, 

ii) (3 = a - the petals degenerate to segments (Fig. E] (b)), 

iii) P G (a, | — at) - the solution does not exist because the map is not conformal; 
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1.0 
0.8 
0.6 
0.4 
0. 




-0.5 0.0 

(a) a = f,/3 



0.5 

_7T_ 

16 




-0.5 0.0 

(c) a=f, = 



0.5 
16 



1 .4 




1 2 




1 .0 




0.8 


1 


0.6 













-1.0 



-1.0 



(b) a = f, (3 

1.4 r 




(d) a = ff , /3 = § 



Figure 5: Two-petal solutions 



B) a = |, /3 G (0, |) - the function /(ii>) can be expressed through elementary func- 
tions: 



/(to) = (w + ur 1 ) [ 1 



1/4 



COS 



9^ ^ • 

z arcsm 

U vr / 



(42) 



(to + w 1 ) , 

a typical pattern is shown in (Fig. 0(c)), the solution degenerates at (3 — f . 
C) aG(f,f) 

i) (3 G (a, |) - the petals are shown in Fig. [5] (d), 

ii) (3 = | — a - the petals degenerate to segments, 

iii) (3 G (| — a, a) - the solution does not exist because the map is not conformal. 

Let us also note that the shape of two-petal solutions at small (3 becomes close to one- 
petal solutions with the same a (Fig. [6]), as it could be expected from the differential 
equations. 
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-0.5 0.0 0.5 

Figure 6: One-petal solution at a — f and two- 
petal solution at a — | and f3 = j^. 
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